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摘 要 : 针对 传统 的 莱 斯 K 因子 一 二 阶 矩 估计 方法 因 贝 塞 尔 鸥 数 存在 计算 复杂 度 大 、 实 用 性 低 的 问题 ， 提 出 一 种 基于 
贝 塞 尔 函 数 阶 数 的 菜 斯 因子 答 估 计算 法 。 该 算法 首先 根据 原 矩 估计 公式 计算 开 值 对 应 的 天 估计 值 ， 然 后 根据 不 同 阶 数 
的 贝 塞 尔 近似 公式 计算 矩 估 计 值 并 求解 对 应 的 开 值 ， 计 算 估计 准确 率 ， 进 而 确定 贝 塞 尔 也 数 阶 数 的 选择 ， 简 化 矩 估 计 
算法 。 实 验 研究 表明 ， s LUN 在 KK 值 较 小 效果 更 优 。 
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Improvement of moments estimation for Rice factor based on Bessel functions 


He Yigang, Liu Chuhuan, Yuan Lifen 
(School of Electrical Engineering & Automation, Hefei University of Technology, Hefei 230009, China) 


Abstract: Aiming at the problem of high complexity and low practicability caused by the Bessel functions in the traditional 
first- and second-order moments estimation algorithm, this paper proposed an improved moments estimation method based on 
the order selection of Bessel function. Firstly, this algorithm calculated the moments estimation value corresponding to different 
K by using original estimation formula, then calculated moments estimation value with Bessel function of different orders, 
further solved the corresponding K value and its estimation accuracy, finally, determined the order of Bessel function acquired 
and improved the estimation method. Experimental results show that the improved method obtains good performance in time 
consuming, especially for small K. 
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出 基于 时 间 序 列 数据 的 天 DS RB f VEU AE: 文献 [8] 对 文献 [7] 

的 方法 的 有 效 性 进行 验证 ， 并 研究 了 环境 对 天 因子 的 影响 ; 文 
莱 斯 信道 是 无 线 通信 信道 常用 的 一 种 典型 信道 模型 ， 用 来 献 [9] 对 文献 [7] 中 的 方法 深入 研究 , 将 背景 噪声 和 干扰 功率 纳入 
描述 接收 信号 的 统计 时 变 特性 目 。 莱 斯 信道 的 接收 信号 包括 直 “模型 ,得 到 一 个 改进 的 矩 估计 算法 , 提高 了 天 因子 的 估计 精度 ; 
视 路 径 分 量 和 散射 路 径 分 量 让。 莱 斯 因子 K 体现 了 莱 斯 信道 文献 [10] 根 据 文 献 [6] 中 近似 方程 导出 玉 因 子 的 一 个 显 式 的 矩 估 
直 视 分 量 功率 和 散射 分 量 功率 的 比例 关系 ， 是 表征 信道 衰落 程 。 计 方 法 ; 文献 [11~13] 提 出 一 种 二 、 四 阶 矩 的 天 因子 闭合 估计 。 
度 和 链 路 预算 的 重要 参数 B44。 莱 斯 因子 越 大 ， 信 道 质量 越 好 ， ”文献 [11，12] 仿 真 表明 文献 [6] 提 出 的 一 、 二 阶 包 络 矩 估计 方法 
因此 莱 斯 因子 K 的 提取 对 分 析 信道 模型 具有 非常 重要 的 意义 。 “在 归 一 化 偏差 和 渐进 方差 性 能 优 于 二 、 四 阶 矩 估计 方法 ， 文 献 
基于 信号 包 络 的 矩 估计 方法 是 提取 信道 莱 斯 因子 的 重要 方 ” [13] 比 较 这 两 种 抢 估 计 方 法 ， 存 在 下 面 异 同 : 二 、 四 阶 矩 的 估 
法 。 这 类 方法 原理 简单 ， 不 利用 相位 信息 ， 不 需要 莱 斯 信道 的 计 方 法 耗 时 稍 少 于 一 、 二 阶 矩 的 方法 ， 一 、 二 阶 矩 的 估计 方法 
先 验 知识 ， 大 大 降低 了 计算 复杂 度 ， 对 算法 的 健壮 性 具有 很 大 ”精度 要 高 于 二 、 四 阶 矩 的 估计 方法 。 其 主要 原因 在 于 : 一 、 二 
t 势 。 近 年 来 ， 很 多 研究 人 员 都 在 天 因子 的 矩 估计 方面 做 了 大 阶 矩 估计 方法 中 含有 贝 塞 尔 函 数 , 函数 本 身 阶 次 高 `. 计 算 复杂 
量 贡 献 和 研究 。 文 献 [5] 提 出 基于 包 络 的 二 、 四 阶 矩 估计 方法 ， — 且 不 能 用 初等 函数 计算 天 因子 的 解析 形式 ， 而 二 、 四 阶 方 法 表 
文献 [6] 提 出 基于 包 络 的 一 、 二 阶 矩 估计 方法 ， 这 两 种 方法 通过 —— 达 式 整体 阶 次 较 低 ， 能 够 求解 出 天 因子 的 解析 形式 ， 导 致 二 、 
将 天 的 函数 等 价 于 所 测 的 包 络 矩 的 比率 来 求 取 天 值 :文献 [7] 提 四 阶 矩 方法 耗 时 小 于 一 、 二 阶 矩 ;其 次 一 、 二 阶 矩 估计 是 基于 
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录用 稿 


对 时 间 序 列 数据 的 一 阶 和 二 阶 矩 计算 , 使 得 莱 斯 包 络 精度 很 高 
因此 ， 如 果 能 够 得 到 一 种 高 效 的 贝 塞 尔 函 数 拟 合 方法 ， 则 能 在 
不 影响 一 、 二 阶 矩 估计 方法 的 计算 精度 的 前 提 下 降低 耗 时 ， 提 
高 算法 的 效率 .因此 , 本 文通 过 对 贝 塞 尔 函 数 进行 多 项 式 拟 合 
提出 了 改进 的 基于  OMER EUIS UN DER c EA 
取 方 法 。 利 用 莱 斯 分 布 信号 的 一 二 阶 和 矩 与 莱 斯 及 因子 之 间 的 函 
数 关 系 ， 结 合 无 线 信道 莱 斯 玉 因 子 的 取 值 范围 先 验 值 ， 以 此 通 
近 贝 塞 尔 函 数 ， 降 低 贝 塞 尔 函 数 阶 数 ， 进 而 降低 整个 表达 式 的 
阶 数 ， 提 高 算法 效率 。 


1 tA A 


131 莱 斯 信道 
莱 斯 信道 中 直 视 路 径 的 信道 系数 并 不 随时 间 变 化 ， 在 信号 
入 射 波 满足 各 角度 均匀 入 射 的 条 件 下 ， 其 概率 密度 函数 为 


| Cor umm ii 


f, o (0.0) RE 


20° 


其 中 : 7,0,@ 分 别 表 示 信 道内 信号 的 衰落 幅度 、 各 径 的 相位 和 


直 视 路 径 的 相位 ，A 为 视 距 路 径 分 量 的 包 络 ， 利 用 上 式 对 相位 
进行 积分 可 得 到 信道 增益 的 包 络 分 布 函数 : 
A? A 
pd er) -— exp( H ED 20 Q) 


人 数 ， 莱 斯 分 布 随机 变量 
(0397; 467 ,? — 42 +2c2，42、2o2? 分 别 表示 了 视 距 分 量 功率 
和 漫 射 分 量 功率 。 信 道 的 莱 斯 因子 定义 为 上 

K = A? /20° Q) 


则 式 (25 可 以 改写 为 
2 

pdf r) = — expl -C +K) 2Kr/o),r20 (D 
o o 


1.2， 莱 斯 因子 的 矩 估 计 法 
莱 斯 分 布 变 量 的 各 阶 矩 表示 如 下 ™: 

HA, = E[r']-7 (o^)? T(n/ 241) exp( CK), F (n/ 2 Ll; K) (5) 
Ep: DO 是 伽 玛 函数 ，, 忆 为 合流 超 几何 函数 。 可 以 看 出 矩 估 
i52) K f oct Xs 因此 要 确定 玉 的 值 , 矩 估 计 的 方法 要 求 
至 少 估计 出 地 的 两 个 不 同 的 矩 ， 常 用 的 有 一 阶 矩 、 二 阶 矩 和 四 
WB. HX CO 计算 可 得 

u -Er]-z(o^"exX-K)FQG/2LK) (6 


[mw 


-Er]-o^-A =0°(K +1) O) 
= E[r^] 2 20^ - 40? A? + A* (8) 
利用 式 (6) (7) 可 得 一 、 二 阶 矩 与 玉 的 关系 式 ; 
-K K.» 
fK) = I: mg DIC KOP © 


EA, CME RAER 1 所 示 。 
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开始 


提取 已 存储 的 信号 
序列 ， 求 其 包 络 


计算 包 络 信号 的 一 阶 敌 / 1、 
Z BME > 


l 


利用 公式 计算 出 K 不 同时 对 应 的 20 值 记 为 A， | 


建立 A 值 与 K 值 对 应 的 表格 


l 


计算 fio (K) — ap ans, NE 


近 的 A 值 ， 则 A 值 对 应 的 K 值 即 为 估计 结果 


图 1 一 二 阶 和 矩 估 计算 法 流程 


2 ” 贝 塞 尔 函 数 多 项 式 拟 合 阶 数 的 选择 


21 贝 塞 尔 函 数 近似 公式 变量 的 取 值 范围 

ER (9) 中 ,由 于 含有 贝 塞 尔 函 数 ， 不 能 用 初等 函数 计算 
玉 因 子 的 解析 形式 ， 使 得 一 二 阶 矩 估计 的 计算 量 大 大 增加 ， 时 
效 性 差 ， 因 此 在 实际 应 用 中 较 少 使 用 ， 无 法 发 挥 一、 二 阶 矩 的 
估计 方法 精度 优势 。 若 对 贝 塞 尔 函数 进行 多 项 式 拟 合 ， 在 精度 
损失 很 小 时 ， 则 能 够 有 效 降低 一 、 二 阶 矩 估计 的 复杂 度 ， 而 损 
失 的 精度 不 影响 估计 性 能 。 


第 一 类 修正 贝 塞 尔 函数 表达 式 为 
x?" x? xí 
Iy(x =]+—+— + 
(x)= D D? 4 64 
x?" 3 5 


X f X ; X 
42"WO-2) 2 


1007 OY. e 


其 中 : 工 为 伽 玛 函 数 。 
文献 [15，16] 提 出 如 下 公式 近似 7 RI T, Cx) : 


2 


4 
1+ 二 + 二 ,0<x<2.5 
Ix) 4 64 


e" [19x,x » 2.5 
x 
一 十 一 十 一 一 ,0 过 x<2.5 


2 16 384 
LOG) . 


X> 2.9 


e 
2.5 x 
Iy Go 的 近似 式 和 理论 值 进行 了 比较 , fO € x € 2.5 d 
习 内 时 ,近似值 与 理论 值 几乎 重合 , 误差 可 以 忽略 ，T1(x) 的 近 
以 式 和 理论 值 比 较 ， 同 样 在 0<xz<2.5 时 ,两 者 几乎 重合 。 因 
Jb, 如 果 变 量 研究 范围 为 0<x< 2.5 , 则 当 贝 塞 尔 近似 函数 的 
阶 次 选取 合适 的 情况 下 ， 可 以 用 式 C110 和 《12) 代替 贝 塞 尔 
函数 对 莱 斯 K 因子 进行 估计 。 
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录用 入 
x? x* 
I(x)214-——-—-.,0€xx2.5 (11) 
4 64 
3. 
LG)-IÍ.2-4..,0xxx2.5 (12) 
2 16 


将 式 (11) C120. 代入 到 式 9) +, HIERBA K 的 范 


和 是 OQ<K<5， 因 此 ， 
0<K<5 范 围 内 时 , 如 果 阶 次 选取 合适 


, 则 完全 


当 实 际 莱 斯 信道 的 莱 斯 


KK 因子 在 


可 以 利用 贝 塞 


尔 函 数 的 近似 函数 式 取代 贝 塞 尔 函 
提取 。 

2.2” 贝 塞 尔 函 数 近似 公式 阶 数 的 选择 
贝 塞 尔 近 似 函 数 的 阶 次 对 估计 精度 存在 影响 


cm 


数 , 简化 莱 斯 信道 K 因子 的 


通常 阶 次 


越 高 ， 精 度 越 高 ， 算 法 也 越 复杂 。 为 保证 莱 斯 k 
足 信道 特征 提取 性 能 需求 ， 并 尽 可 能 的 降低 算法 
何 选 出 最 佳 贝 塞 尔 近似 函数 的 阶 次 显得 十 分 重要 
的 贝 塞 尔 近似 函数 阶 次 的 选择 流程 如 图 2 所 示 。 


AT EUH ES 
的 复杂 度 ， 如 
。 据 此 ， 本 文 
具体 表述 为 : 


a) 计算 基于 贝 塞 尔 函 数 的 不 同 k 
估计 值 ; 

b) 计算 所 选取 阶 次 的 贝 塞 尔 近 
计算 对 应 的 莱 斯 因子 估 R 值 ; 


似 函 数 一 


值 时 所 对 应 的 一 、 


、 二 阶 矩 ， 


二 阶 矩 


并 以 此 


c) 比较 理论 天 值 与 估 R 值 ， 计 算 莫 


相对 误差 ， 如果 误差 在 


es ds 
加 1， 然后 重复 步骤 b 

图 3 和 4 分 别 是 零 E 阶 贝 塞 尔 
与 理论 值 的 比较 , 由 图 3 TARH K ELO, 5] 内 ， 
差 较 大 , 4 阶 和 6 阶 能 较 好 跟踪 理论 值 , 图 4 中 ， 
理论 值 较 小 ，5 阶 和 7 阶 基本 上 与 理论 值 重合 。 


GD 


Y 


Te X 


函数 近似 式 
2 阶 近似 式 误 
3 阶 近 似 式 与 


2N EKI 


利用 式 (8) HAHKETENA EN- MEHA fiK) 


记 为 A 


v 


NEREAREN I(x) a T, (x) 阶 数 的 选取 n1, n2 


v 


n1=n1+1;n2=n2+1 


Y 


利用 贝 塞 尔 函 数 近似 式 计算 K 值 不 同时 对 应 一 


二 阶 矩 估计 值 ， 记 为 B 


v 


建立 A 值 和 B 值 对 应 的 K 值 的 表格 


y 


— 10(x) 
rus Xx<2.5 的 2 阶 近 似 
一 一 x<2.5 的 4 阶 近似 
2 * x<2.5 的 6 阶 近 似 


Bessel 函 数值 


图 3 不 同 零 阶 贝 塞 尔 函数 近似 式 与 理论 值 曲线 
3 
25. | 一 1100 
-一 X<2.5 的 3 阶 近 似 
* X<2.5 的 5 阶 近 似 
im ?| ”|-5-x<2.5 的 7 阶 近似 
| 
2.5 
图 4 不 同一 阶 贝 塞 尔 函 数 近似 式 与 理论 值 曲 线 


3 ”实验 仿真 分 析 


K 有 效 性 评价 
为 了 对 贝 塞 尔 函 
有 效 性 评价 ， 
的 有 效 性 进行 评价 : 


3.1 


a) 1 


EHZ D: D= 7 x100% 


的 取 值 长 度 ，m 表示 选取 阶 数 分 别 nl,n2 
E K KEX L 内 估计 的 包 值 与 原 一 二 阶 和 
相同 的 个 数 。 


N 


D) SEZTRERSZS RMSE gp 


i=1 


数 近 似 式 选取 不 同 阶 数 计算 得 到 不 值 进行 
本 文 主要 以 准确 率 作为 对 下 的 不 同 取 值 长 度 估算 


BENE E 


贝 塞 尔 函 


数 近 似 公 


斯 因子 天 


N 式 


EE 估计 得 全 


` (X, E Au H 


| K 值 完全 


中 和 


de BUB KE L 内 的 样 
表示 第 i 个 理论 值 。 
3.2 仿真 结果 

本 文 第 一 类 零 阶 修正 贝 塞 尔 函 
类 一 阶 修正 贝 塞 尔 函 


AUR, XQ 表示 第 i 


N 


个 样 


本 观测 值 , X mi 


数 阶 数 选择 为 4 和 6, , 第 一 


数 阶 数 选择 3 和 5， 对 其 进行 组 合 ， 得 到 


从 B 中 查 到 与 A 中 最 接近 的 值 对 应 的 K 值 ,计算 K 值 在 长 度 内 的 准确 率 


图 2 贝 塞 尔 近 似 公式 阶 数 选取 流程 


第 一 类 零 阶 修正 贝 塞 尔 函 数 阶 数 为 4、 


第 一 类 一 阶 修正 贝 塞 尔 


函数 阶 数 为 3， 


第 一 


KEH 


阶 修正 贝 塞 尔 函 数 阶 数 为 5， 


FE 贝 塞 尔 函 数 阶 数 为 4、 第 一 类 一 
第 一 类 零 阶 修 正 贝 塞 尔 函数 阶 数 


为 6、 第 一 类 一 阶 修正 贝 塞 尔 函 
贝 塞 尔 函数 阶 数 为 6、 第 一 类 一 阶 修正 贝 塞 尔 函 
4 种 组 合 ， 具 体 公式 如 下 : 


数 阶 数 为 3 和 第 一 类 零 阶 修正 
数 阶 数 为 5 共 
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3048 ITER, F: AFRA AN AFE kt 
i x x 根据 求 得 的 近似 矩 估 计 值 求解 对 应 的 莱 斯 因子 玉 值 ,与 原 
I =]+— +— ; 
组 合 1 YQ) 4 64 0<x<25 K (RECS IEF Kog Kiss Kies Kiss 分 别 是 组 合 1、2、3、4 对 
TR TEE din 应 估计 得 到 的 KK 值 , 表 1 给 出 了 四 种 组 合 天 不 同 取 值 范围 的 均 
X) 三 一 十 一 
: 2 16 方 根 误差 。 
xX ox 表 1 4 组 合 对 应 的 不 同 取 值 范围 的 均 方 根 误差 
1,(X)=1+—+— K K K K 
组 合 2 4 64 0<x<2.5 1243 1245 1263 1265 
x x x (0,2) 0 0 0 0 
I(x) = 一 + 一 十 一 一 
2 16 384 (0,4) 0.5760 0.0643 0.3539 0 
x o (0,5) 2.7760 | 0.6975 1.8339 0.0231 
ROO oor 0301 表 2 4 种 组 合 对 应 的 不 同 取 值 范围 的 准确 率 (%) 
组 合 3 ! ,0<x<2.5 ! 
L (x) E K Ks Kys Kos Kiyes 
2 16 (0,2) 100 100 100 100 
"AE. (0,4) 52.38 66.67 52.38 100 
Boc rt (0,5) 4231 53.85 4231 84.62 
EN 4 64 2904. «25 ， : : : i 
Peta c 由 图 6、 表 1 和 2 可 以 看 出 , 在 0<KK<2 区 间 内 , 组 合 1、 
(x) = 
2 16 384 2、3、4 计算 得 到 的 玉 均 方 根 误差 为 0， 准确 率 均 为 100%， 在 


分 别 对 组 合 174 得 到 的 贝 塞 尔 近似 公式 计算 出 不 同 值 时 对 0< 玉 <4 区 间 内 , 组 合 1 均 方 根 误差 值 最 大 , 组 合 3 均 方 根 根 
应 一 二 阶 矩 估计 值 以 及 原 一 二 阶 矩 公式 计算 出 不 同 值 时 对 应 的 误差 次 之 ， 这 两 种 方法 的 准确 率 相 同 ， 组 合 2 估计 得 到 天 值 的 
一 二 阶 矩 值 ,其 中 fi, 表 示 原 始 一 二 阶 矩 估计 公式 对 应 不 同 玉 。 均 方 根 误差 和 准确 率 优 于 组 合 1 和 3， 组 合 4 均 方 根 误差 为 0， 
时 计算 值 ，fi2s 表 示 组 合 1 对 应 不 同时 的 计算 值 ，fiwss 表 示 。 估计 准确 率 为 100%, 在 该 区 间 内 能 够 完全 准确 估计 出 天 值 , 在 
组 合 2 对 应 不 同 久 时 的 计算 值 ，fizes 表 示 组 合 3 对 应 不 同时 。 区 间 0< K «5 A, 相 比 前 两 个 区 间 , 均 方 根 误差 都 相对 增 大 ， 
KHE fis ERES 4 对 应 不 同 尺 时 的 计算 值 ， 仿 真 ”准确 率 相对 下 降 ， 组 合 4 优 于 其 他 组 合 1、2、3。 
KK 值 的 步 长 取 0.2， 仿 真 结果 如 图 5 所 示 。 33 ”复杂 度 分 析 
图 5 可 以 看 出 ， 不 同 组 合 的 值 与 原 一 二 阶 计算 值 的 精度 组 合 1 4 的 最 高 阶 数 分 别 为 4、5、6、6， 本 文 为 保证 贝 守 
随 着 K 取 值 范 围 长 度 变 化 而 改变 ,在 区 间 0< < 2 内, 组 合 — 尔 函 数 的 计算 精度 ， 选 取 的 阶 数 为 1022， 天 步 长 为 0.2, 在 
1、2、3、4 计 算 值 与 原 一 二 阶 计算 值 基本 相等 ,在 区 间 0<K<4 0<K<2 区 间 内 , 组 合 1 4 计算 得 到 的 均 方 根 误差 为 0, 组 合 
内 ， 组 合 4 计算 值 与 原 一 二 阶 矩 计算 值 一 致 ， 组 合 2、3、1 计 1 的 阶 数 最 低 ， 在 估计 效果 相同 的 情况 下 选择 组 合 1 来 近似 一 
算 精度 依次 降低 ， 不 能 满足 精度 要 求 。 随 着 天 取 值 的 增 大 , 组 ”二 阶 矩 估计 公式 ， 最 高 阶 数 为 4， 在 很 大 程度 上 降低 计算 复杂 
frl. 2. 3. 4 曲线 偏离 理论 值 ， 组合 1 偏离 最 大 ,组合 4 的 偏 。 度 ; 在 0< 玉 <4 区 间 内 , 组 合 1 和 3 两 种 方法 的 准确 率 相同 ， 


BUR Ait. 但 组 合 1 均 方 根 误差 较 大 ， 组 合 2 估计 得 到 均 方 根 误差 和 准确 
- 率 优 于 组 合 1 和 3， 天 值 稍 大 于 2 时 可 以 选择 使 用 组 合 4， 组 
| 合 4 均 方 根 误差 为 0， 估计 准确 率 为 100%， 在 该 区 间 内 能 够 完 


全 准确 估计 出 值 , 选择 组 合 4 作为 该 区 间 的 一 二 阶 矩 估计 公式 ， 
阶 数 为 6 阶 ， 降 低 计算 复杂 度 。4 种 组 合算 法 及 原始 算法 耗 时 


如 表 3 所 示 。 
表 3 4 种 组 合算 法 及 原始 算法 耗 时 (ms) 

K K, Ks Kiys K263 Koss 
(0,2) 4.0 0.0794 0.0862 0.0929 0.1013 
(0,4) 4.2 0.0964 0.1014 0.1087 0.1212 
(0,5) 4.3 0.1023 0.1073 0.1228 0.1269 


表 3 可 以 看 出 ，4 种 组 合算 法 的 计算 时 间 明 显要 小 于 原 
一 二 阶 矩 估计 算法 运行 时 间 ， 根据 KK 区 间 选 取 贝 塞 尔 函 数 的 阶 
数 ， 降 低 矩 估计 算法 的 计算 复杂 度 ， 实 现 算法 的 优化 改进 。 


4 ”结束 语 
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6 4 种 组 合 求解 出 的 莱 斯 天 值 


针对 一 二 阶 矩 估计 公式 中 贝 塞 尔 函 数 阶 次 较 高 导致 计算 复 


录用 稿 


杂 度 高 的 问题 ， 提 出 一 种 基于 贝 塞 尔 函 数 近似 式 阶 数 的 选取 方 
法 的 莱 斯 天 因子 的 一 二 阶 矩 估计 算法 。 通 过 对 本 文 构建 的 4 种 
不 同 组 合 仿真 实验 ， 验 证 了 天 值 取 值 区 间 在 0-4 内 能 够 保证 KK 
估计 准确 的 前 提 下 贝 塞 尔 近 似 式 阶 数 选取 小 于 理论 选择 的 阶 数 ， 
减少 了 算法 的 运行 时 间 ， 提 高 了 算法 时 效 性 ， 实 现 了 算法 的 改 


t. 
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